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G. BARLES(i), F. DA LIO^^), P.-L. LIONS^^) AND P.E. SOUGANIDIS^^) 

Abstract. We study periodic homogenization problems for second-order pde in half- 
space type domains with Neumann boundary conditions. In particular, we are interested 
in "singular problems" for which it is necessary to determine both the homogenized 
equation and boundary conditions. We provide new results for fully nonlinear equations 
and boundary conditions. Our results extend previous work of Tanaka in the linear, 
periodic setting in half-spaces parallel to the axes of the periodicity, and of Arisawa in 
a rather restrictive nonlinear periodic framework. The key step in our analysis is the 
study of associated ergodic problems in domains with similar structure. 



1. Introduction 

We study issues related to the homogenization and ergodic problems for fully nonlinear, 
non-divergence form, elliptic and parabolic boundary value problems in half-space type 
domains with possibly nonlinear Neumann boundary conditions. In particular, we are 
interested in problems for which it is necessary to identify both the homogenized equation 
and boundary conditions. Our results represent a first step towards the resolution of the 
problem in general domains, which remains open even for the linear non-divergence form 
problem. The situation is, of course, different for divergence form equations, since, in that 
case the boundary conditions are encoded in the variational formulation. 

In order to be more specific and to describe the problem in a clear way, we first discuss 
heuristically two model cases involving linear elliptic equations which lead to different 
behaviors and difficulties. To simplify the presentation, we assume that all the functions 
appearing in the examples below have the needed regularity properties and are periodic 
with respect to the fast variable x/e, which we denote throughout the paper by y. 
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The first problem, which we call "regular" , is 

. ^ J —ti{A{x,e^^x)D'^u^) — b{x,e~^x) ■ Du'^ + u'^ = f{x,e~^x) in (9 , 

^ > \ Du'~j = g on dO, 

where O is a domain in M.^ , not necessarily of half-space type, Du"^ and D'^u'^ denote 
respectively the gradient and Hessian matrix of the solution m^, and 7 and g depend only 
on x. 

The typical "singular" problem is 

. ^ J —ti{A{x,e^^x)D'^u'^)—E^^b{x,e~^x)-Du^ + u'^ = f{x,e~^x) in O , 
^ \ Du' -J = g on dO , 

which has as a particular case the divergence form equation 

-div(A(x, e~'^x)Du^) + = f{x, e~^x) in O . 

If we use the formal expansion 

ti^(x) = u{x) + ev{x, e~^x) + e'^w{x, e~^x) + 0{e^) , 

then, in the "regular" case, we find that the leading (the e~^) term yields 

-tT{A{x,y)Dlv) = inR^; 

a standard Liouville-type property implies that v is independent of y. Moreover, when 
7 and g depend only on x and not on y, the function u is expected to satisfy the same 
boundary condition as m^, since the formal expansion yields 

Du" ■ 7 = ■ 7 + eD^v ■ 7 + eDyW ■ 7 + 0{e'^) . 

This problem can be treated, at least for uniformly elliptic A, using the perturbed test- 
function method of Evans p2]- If, however, either 'y or g depend on y, different arguments 
are needed. 

The same expansion, in the "singular" problem leads ot 

tT{A{x, y)DlyV) + b{x, y) ■ {Du + Dyv) = , 

where now the function v actually depends on y in general. This interferes with the 
boundary condition, since the expansion is now 

Du^ ■ 7 = ■ 7 + DyV ■ 7 + 0{e) , 

and clearly v plays a role in determining the boundary condition for the limiting equation. 

This is the main problem in the singular case. The key issues are the identification of 
the equation and the boundary condition for u. In the "regular" case, the issue is the 
study of the asymptotic limit for y-dependent 7 and g. 

There exists an extensive body of work dealing with the homogenization of the Dirichlet 
problem for fully nonlinear first- and second-order partial differential equations in periodic, 
quasi-periodic, almost periodic and, more recently, stationary ergodic media. Listing 
references is beyond the scope of this paper. 
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Little is, however, known for the Neumann problem except for divergence form equa- 
tions with the usual (co-normal) boundary conditions, which are treated in the classical 
book of Bensoussan, Lions and Papanicolaou [9j. 

In [19], Tanaka considered the two model problems discussed earlier by purely proba- 
bilistic methods in half-spaces. In pQ, Arisawa studied special cases of homogenization 
problems again in half-space type domains under rather restrictive assumptions. Our 
methods are inspired from [T] but yield more general results. 

To study the asymptotic behavior of the m^'s, we first consider the usual cell (ergodic) 
problem which is supposed to give the equation inside O. In the "singular" case, it is 
formulated as follows: 

For each x,p & R^, find a unique constant X{p,x) such that there exists a bounded 
solution V of the equation 

(1.3) -tT{A{x,y)DlyV)-b{x,y)-{p + Dyv) = X{p,x) in . 

The map p i-^ X{p,x) is clearly linear. Moreover, the only interesting case is when 
X{p,x) = 0, otherwise we get a trivial first-order equation for u. We do not discuss here 
the type of conditions on A and b which yield X{p,x) = 0. We just point out, however, 
that, for divergence form equations, this is true, as it can be easily seen by integrating 
the equation over a period. 

For the boundary behavior, we restrict, for simplicity, to the half-space case. To fix the 
notation, we assume that O = {xn > 0}, in which case O = e^^O. The fact that both 
O and do are invariant under the scaling is a property which is needed, at least for our 
approach, in the study of general domains. 

The natural ergodic problem on the boundary is to find, for each p,x E M^, a unique 
constant Ji{x, p) such that there exists a bounded solution w of the boundary value problem 




-tr(A(x, y)DlyW) - b{x, y)-{p + Dyw) = in O , 
{p + Dyw) ■ 7(0;, y) = g{x, y) + /i(x, p) on dO . 



This type of problem was first studied in [T] both in bounded and in half-space type do- 
mains for Hamilton- Jacobi-Bellman type equations with oblique Neumann boundary con- 
ditions but under rather restrictive assumptions. We provide here existence and unique- 
ness results for general nonlinear equations and some boundary conditions. 

The paper is organized as follows: In Section 2 we introduce the hypotheses and state 
and prove the existence result for the nonlinear boundary ergodic problem. Section 3 is 
devoted to the study of the uniqueness properties of /i(x,p). In Section 4, we consider 
homogenization problems. We conclude with some remarks about the non-periodic case. 



2. Existence results for ergodic boundary problems in half-spaces 



We study the following boundary ergodic problem : Find a constant fj, such that there 
exists a bounded, continuous solution of 



(2.1) 



F(D^u,Du,x) = A in O, 
L{Du,x) — IjL on do. 



Here O C is a smooth, periodic, half-space-type domain (see (01) below for the 
definition), F and L are real valued, at least, continuous functions defined on S'^ x x O 
and x dO respectively, where is the space oi N x N symmetric matrices, and A is 
a constant which is chosen in a suitable way below. 

We formulate next the assumptions on F, L and O. Note that, for the sake of the clarity 
of the exposition, in this note, we do not attempt to state the most general conditions. 

As far as F is concerned, we assume that: 

(FO) F : xR^ xR^ ^Ris Si continuous, Z^-periodic function in x, 

F is locally Lipschitz continuous and there exists a constant K > 



(Fl) 



such that, for all x,y eR^ , p,q eR^ , and M,N e , 

\F{M,p,x) - F{N,q,y)\ < K{\x - y|(l + \p\ + |g| + ||M|| + ||iV||) 



+ \p-q\ + \\M-N\\} 



(F2) 



, p e 



\{0}, 



there exists k > such that, for all x e 
and M,N eS^ with N > 0, 

F{M + N,p,x) -F{M,p,x) < -K{Np,p) + o{l)\\N\ \ , 
with o(l) — > as |p| — > oo, 
where, for peR^\ {0}, p = p/\p\. 

As far as the domain is concerned, we assume that 

(01) C is a smooth, periodic half-space type domain, 

i.e., that: 

(i) do is a C^-boundary, 

(ii) there exists {N — 1) hnearly independent vectors /i, • • • /at-i £ '^^ , such that, for all 
X e dO, X E O and ^i, • • • -Zat-i £ Z, 

X + zifi H h ZN-ifN-i e dO and x + Zifi H h z^-ifw-i e O , 

and 

(iii) there exists a unit vector /jv G R^, orthogonal to /i, • • • , /at-i, and ^ > such that, 
for all X e O and R> R, 

X + R/n e O . 



Examples of domains satisfying (01) are 

O = {x = {Xi, ■ ■ ■ ,Xn) '■ Xat > 0} 

or, more generally, 

O = {x = (xi, • ■ ■ , Xn) : Xn > ^(xi, ■ ■ ■ , Xn-i)}, 
where ifj G ^^(M^-i) is Z^-^-periodic . 

We denote by d the sign-distance function to dO, normalized to be positive in O and 
negative in \ O, and we recall that, for all x G dO, the outward normal n{x) to dO 
at X is given by 

—n{x) = Dd{x) . 

A key ingredient in the existence proof are the up to the boundary C^'^-regularity 
results for Neumann boundary value problems which were obtained by Barles and Da Lio 
in [5]. To use them, it is necessary to introduce the following set of assumptions on L. 

There exists u > such that, for all {x,p) G dO x and t > 0, 
L{p + tn(x), x) — x) > ut . 



(LI) 



(L2) 



There exists constant K > such that, for all x, ?/ G dO and g G 
\L{p,x)-L{q,y)\<K[{l + \p\ + \q\)\x-y\ + \p-q\] . 



There exists a G C(]R^ x O) such that, 

(L3) ^ as t ^ +00 and locally uniformly in {p,x), 

t-^L{tp,x) L^{p,x) . 

The regularity results of [5] depend on whether L is linear or nonlinear and require 
some additional assumptions on F. In the nonlinear case, it is necessary to strengthen 
(F2) and to ask for F to be a uniformly elliptic, i.e., to assume that 

there exists k > such that, for all x G O, p G M^, and M,N e S'^ with > 0, 

F{M + N, p, x) - F{M, p, x) < -Kti{N) . 

Moreover, we need to assume that 

there exists F^o G C{S'^ x x O) such that, locally uniformly in (M,p, x), 

ast->oo, t-''F{tM,tp,x) ^ F^{M,p,x) 

For linear boundary condition, i.e., if 

L{p,x) = p-7(x) - g{x) , 

with 7 G C°'i(9C; M^) and g G C°'^((90) for some (3 G (0, 1), instead of (F3) we assume 
that 



(F3) 



(F4) 
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(02) I ^^^^^^ ^ ^ C'°'^(C; S'(A^)) and cq > 0, such that 

A > cqIcL , and yl(x)7(x) = n{x) on (90 , 
and, for all i? > 0, 

{there exist Lr, A^j > such that, for all x G O, \u\ < R, \p\ > Lr 
and M,M E with M > 0, 
F(M + M,p,x) - F{M,p,x) < -\r{MA^^^)p, A^^^)p) + o{l)\\M\\ , 

where o(l) ^ as \p\ oo. 

We call (F5) "adapted ellipticity" , since the ellipticity condition is adapted in the 
direction of the oblique derivative through A. 

We remark that all the assumptions on F are satisfied by Hamilton- Jacobi-Bellman or 
Isaacs Equations under the usual conditions on the coefficients, while (LI), (L2), (L3) are 
natural conditions for (nonlinear) Neumann boundary conditions. 

The constant A in fl2.ip is given by the next proposition. Since the proof follows the 
argument of Barles and Souganidis [S], we omit it. 

Proposition 2.1. Assume (FO), (Fl) and (F2). There exists a unique constant A such 
that there exists a periodic solution u G C°'^(]R^) of 

(2.2) F{D%Du,x) = \ mR^ . 

Our result is 

Theorem 2.1. Assume (01), (FO) and either (Fl), (F3), (F4) and (LI), (L2), (L3) if L 
IS nonlinear or (Fl), (F2), (F5), (02), 7 G C^'^{dO) and g G C^'^{dO), if L is linear. 
There exists fi such that (12.1 p has a continuous bounded viscosity solution, which has the 
same periodicity property as the domain O. 



Proof. We concentrate first on the nonlinear case. To simplify the exposition, we assume 
that /tv = ctv and span(/i,-- - , fN-i) = span(ei,-- - ,eAr_i), where (ei,-- - , cat) is the 
standard orthonormal basis of M^. In addition, we assume that G dO. 
For 0<£:<a<l, we introduce the approximate problem 

. . j F{D'^u, Du, x) + eu = X + eu in O, 

^ \ L{Du,x) + au = on dO, 

where F and u are given by Proposition 2.1. 

The existence and uniqueness of u follows from classical arguments from the theory of 
viscosity solutions, and in particular, the Perron's method [15] and the comparison results 
[3]. Moreover, in view of its uniqueness, the solution has the same periodicity properties 
as the domain, i.e., for all x G O and ki - ■ ■ k^_i G Z, 

u{x + kici H h kN-ieN-i) = u{x) . 



Finally standard comparison arguments (see |3]) yield that 

niax('u — m) < sup |L(e,x)|, 

O xGR^,|e|<||DS||oo 

which, of course, implies that 

au is uniformly bounded in x, a and e ; 

notice that the e-term in f l2.3p is introduced just in order to prove this estimate. 
For i? > sufficiently large, we consider next the domain 

OR = {xeO: Xn <R} 

and solve the new problem 

{F{D^ur, Dur, x) + eur = X + eu in Or, 
Dur ■ = on {xn = R}, 
L{Dur, x) + auR = on do, 

where hr is the external normal to Or on {xn = R}. 

Again it is a standard fact that (12. 4p admits a unique continuous viscosity solution ur 
with the same periodicity properties as the domain. Moreover, the maximum principle 
yields that, for all x E O, 

\ur{x) -ttR(0)| < \\UR-UR{0)\\L-O(^do) , 

and since dO is periodic, we may assume that all the ur — ur{0)^s attain their maximum 
and minimum values at points which remain in a compact subset of M^. 

Letting R oo and using the standard fact that, a.s R ^ oo, ur ^ u in C{0), we 
find that, for all x G (9, 

\u{x) - n(0)| < ||u - m(0)||l°o(oo) • 

We claim next that 

remains bounded as first e — and then a 0. To this end, we argue by contradiction 
and introduce the function w : O ^ M. given by 

u-u(0) 
w = — , . 

I |u-m(0) 11^^(0) 

It is immediate that \w\ < 1, w{0) = and w is a solution of a nonlinear Neumann type 
problem with nonlinear it ies satisfying assumptions (Fl), (F4), (F5), (LI), (L2) and (L3) 
with uniform constants. The C*^'"-estimates of [5] yield that the w's are locally uniformly 
bounded in C°'". Extracting a subsequence, we may assume that, as e — > and a — *• 0, 
the w's converge to w, which solves 

. Foc{D^w,Dw,x) = in C, 

Loo{Dw,x) = on do, 
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and, moreover, 

w(0) = 0, and ||w||oo = l 

with the sup-norm achieved at some point of dO, a fact which contradicts the strong 
maximum principle of [4]. 

Since, by the previous step, u — u{0) remains bounded, we can use again the C°'"-local 
estimates of [5] and pass to the hmit, letting first e — > and then a — > 0. It follows that, 
up to subsequences, 

—au{0) — > /i . 

The proof of the linear case follows the same arguments using a different C'^'^-local 
estimate in [5]. It is worth pointing out that the strong maximum principle of [Ij is still 
vahd under assumption (F5). □ 

3. Uniqueness results for boundary ergodic problems for 
Hamilton- Jacobi-Bellman type equations in half-space type domains 

We present here two proofs for the uniqueness of the boundary ergodic cost fi. The 
first one assumes that F is convex and uses the stochastic control interpretation of the 
problem. It turns out, however, that it is also possible to present a proof which is based 
entirely on pde-type arguments and does not utilize convexity and stochastic control. 

3.1. The stochastic control proof . 
We consider the solution {{Xt)t>o, {kt)t>o) of the stochastic differential equations 



(3.1) 



dXt = b{Xt, at)dt + v^a(Xt, at)dWt - dkt , Xq = x G O, 
h = laaiXs)^{X,)d\kU , XtEO, t>0 



where {Xt)t>o is a continuous process in and {kt)t>o is a process with bounded vari- 
ation. Here {Wt)t>o is an X-dimensional Brownian motion, the control (at)t>o is a pro- 
gressively measurable process with respect to the filtration associated to the Brownian 
motion with values in a compact metric space A and the drift b and the diffusion matrix 
a satisfy the classical assumptions 

b G C(0 xA;R^),ae C(0 x A; S{N)) and, for each a e A, 
(CI) { b{-,a)eC'''\0;R''),ai-,a)eC'''\OxSiN))anda = aa^ 
is uniformly elliptic, with all the constants uniform in a. 

It is known that (CI) yields the existence of a unique solution {{Xt)t>o, {kt)t>o) of 
fl3.ip . (See Lions and Sznitman |18j for the existence and Barles and Lions ^ for the 
uniqueness.) 

The nonlinearity F of the associated Bellman equation and the boundary condition for 
the ergodic boundary problem are given respectively by 

(3.2) F{M, p, x) = sup < — tr [a(x, a)M] — b{x, a) ■ p — f{x, a) 



(El) 



and 

(3.3) L{Du,u,x) = Du ■ 'J — g + fi = ondO. 

We assume that 

{there exist (3 G (0, 1) and u > such that, for all a E A, /(■, a), 
7 and g satisfy : /(-,«) G BUC{R^), g G C°'^(9C), 
7 G ^0.1(90; R^) and 7(x) ■ n{x) ^ on 9C. 
Finally, it is necessary to make the additional "ergodic"-type assumption that 
there exists a bounded, Lipschitz continuous subsolution w of 
F{D^w,x,-fN + Dw,x) = inM^, 
where /at appears in (01) and 

F{X, p, x) = sup I — -tr [a{x, a)X] — b{x, a) ■ p\. 

We remark that, for linear equations, (El) is not a real restriction. Indeed, in the 
Introduction, we explain that X{p,x) = for all p is natural in our framework and the 
condition (El) reads A(— /tv,^) = 0. 

The result is 

Theorem 3.1. Assume (CI), (C2) and (El). There exists at most one constant fi which 
solves the boundary ergodic control problem (12.1 p with A and F and L given by Proposi- 
tion \2.1\ and (13. 2p and (13.31) respectively. 

The key step in the proof of the uniqueness is the following lemma. 

Lemma 3.1. Assume (CI), (C2) and (El). There exists x G O such that 



(3.4) inf 

{at)t 



d\k\ 



+ 00 . 



Proof. Following the proof of Theorem 12. H we assume that = ^n, and we work in 

Or = {xeO : Xn < R} , 

Define 

v{x) = —Xn + w{x) 
where w is given by Proposition 2.1. Then 

F{D^v,Dv,x) = inM^, 

Since v is Lipschitz continuous and w is bounded, the exists a C > such that v 

Df ■ 7 < C on do and v<—R + C on xn = R- 
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Let be the first exit time through the boundary {xat = R} of the process {Xt)t>o 
starting at x. Using the dynamic programming principle for the stopping times we find 
that, for all t > 0, 



v{x) <C inf 

{at)t 



d\k\ 



-iR-C)x{x,t) 



where 



X{x,t) = inf [l{r,<i}] • 

("t)t>0 



Since a is nondegenerate, it follows that, for each fixed R, as t 
uniformly in x, 

which, in turn, implies that, locally uniformly in x, 



oo and locally 



lim inf inf lE^ 

t^oo (at)t>o 



d\k\, 



> C-\R-C + v{x)) 



□ 



Since R is arbitrary, the result now follows for all x G (9. 

We return now to the proof of Theorem 3 which is based on the intuitive idea that, 
to have a unique fi, the boundary needs to be seen in a "sufficient way". This last fact 
which is quantified in a precise way by (13.41) . A similar fact is also a key point for the 
uniqueness of fi in bounded domain. 

Proof of Theorem 13.11 Assume that there exist two constants fii and /X2 such that 
/ii < /i2 with corresponding solutions Ui and U2- The dynamic programming principle or 
the uniqueness for the associated time dependent problem yield, for all t > and i = 1,2, 



UiiX 



inf 

("i)t>0 



[ [f{Xs,a,) + X]ds+ [ [g{Xs)+fii]d\k\, + u,{Xt 
Jo Jo 



Choosing an ^-optimal control for U2 and subtracting the inequalities for ui and U2 we 
find 



ui{x) - U2{x) < {fli - fIi)E^ 

which yields the estimate 



d\k\ 



Uo 



+ UiiXt) - U2iXt) + 6 



(/i2 - fil)E^ 



d\k\ 



Uo 



< 2(||m||oo + ||m||oo) +£ ■ 



□ 



If fi2 > /^i this last inequality contradicts Lemma 3.1. 

It turns out that, for /i to be unique, it is necessary to assume (El). Indeed, consider 
the simple Neumann problem 

—{p" — {p' = in (0, +oo) and — ip'{0) = /i . 

Since (j){x) = //exp(— x) solves this problem for all /i ^ 0, it turns out that since f{x) = 
/iexp(— x), any ^ is a solution of the associated ergodic problem. 
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On the other hand, a direct integration yields that there does not exist a bounded 
subsolution of 

-i)" - V^' + 1 = in (0, +oo) . 

Finally, we remark that in [1] the uniqueness is proved in the uniformly elliptic case 
under the additional assumption that, for all x and a, 

6iv(x, a) < , 

a fact which allows to choose w = in (El). 

3.2. A pde proof and a more general result. 

We state and prove using pde techniques a stronger uniqueness result than the one 
asserted in Theorem 3.1. To this end, we need to assume that there exist 

a homogeneous of degree one, uniformly elliptic, F G C{S^ x x M^) 

(F6) { such that, for all x G M^, p, g G and Mi, M2 G , 

F{Mi,p,x) - F{M2,q,x) < F(Mi - M2,p - g, x) , 

and 

a Lipschitz continuous subsolution v of 
(E2) <( F{D^v,Dv,x) = in O , 

such that, uniformly with respect to x', f (x',XAr) —00 as x^ +00, 

where, for X G M^, we write ) with 

It is immediate that the convex function F given by fl3.2p satisfies (F6), while (E2) is 
a weaker version of (El). We also remark that, although all these assumptions appear a 
bit artificial, [Ij provides some counterexamples. 

The result is 

Theorem 3.2. Assume (F6) and (E2). If F satisfies (Fl), there exists at most one 
constant ^ solving the boundary ergodic control problem. 

Proof of Theorem 13.21 1. Let /ii, fi2, ui, U2 and Or be as in the proof of Theorem 
and assume that /i2 > /^i- Then w = m — m is a bounded supersolution of 

F{D^w,Dw,x) > in O , 
Dw ■ 7 > /i2 — yWi on do. 

For some c > 0, assume that there exists G C{0 x [0, 00)) such that 

(3.5) ^ c on {xjv = M} and 4jr{-, 0) ^ w on dO , 

and 

' + F{D^ip^, Dip^, x) ^ on , 

(3-6) <; 

^— S - on DOr . 
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The comparison principle of viscosity solutions and the fact that w is bounded yield 
the existence of another constant C > such that 

ij^^w + C in C X (0,oo) . 

This last inequality leads to a contradiction, if it is possible to choose the ip^^s so that, 
in addition to (13. 5p and (I3.6I) . they also satisfy, 

(3.7) lim \im ip^(-,t) = +00 locally uniformly in O. 

R^oo t^oo 

For each i? > 0, we define by 

ijj^ = 6v + 6mjiX^ , 
where (5 > is chosen below sufficiently small enough, 

niR = inf {—v{x',xn))^oo as R oo , 

and G C(C X (0, oo)) is such 

(3.8) ^ ^ 1 and lim x^{'^ ^) = 1 locally uniformly in O . 

Such ip^ clearly satisfies (13. 7p and the first inequality in (11). To prove that ip^ is a 
subsolution we compute 

+ F{D'^il)^, D^^, x) = 5mRxf + F(5D\; + 6mRD\^, 6Dv + dniRDx^, x) 

< F{6D^v, 6Dv, x) + 6mR (xf - M^D^^) - C\Dx^\) , 

where is the minimal Pucci's operator associated with F. 

This computation shows that it is enough to choose x^ to be the unique solution of 

'x^-M4D\^)~C\Dxr\=0 in Or x{0,+oo), 

Dx^ -7 = on 9C X (0, +oo) , 

< 

X^ = 1 on {xn = -R} X (0, +oo), 
= on O X {0} , 

which satisfies (13.81) . 

With this choice of x^) ''P^ satisfies the subsolution inequality. The fact that v is 
bounded yields, if 5 is chosen sufficiently small, that the boundary condition in (13.51) and 
the second inequality in (13. 4p are also satisfied. 

□ 

We conclude this section emphasizing that an import issue is to understand the way on 
which /i depends on F and L. This question was addressed in [1] in bounded domains. 
In the context of this paper, the dependence on L is an easy consequence of the existence 
proof, while the dependence on F is not so obvious. 
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4. The homogenization problem in half-space type domains 
with oscillating boundaries 

To make the main ideas of our approach clear, we begin with the hnear problem and 
we consider domains of the form 

Oe^{xeR^ : XN> eilj{e-^x)} = e'^O ^ {y eR^ : vn > ^{y')} , 

where, as before, for x e R-^, x' — (xi, • • • ,xn^i) and ip is a smooth (at least C^'^) 
Z-'^~^-periodic function. 

The model homogenization problem we are interested in is 

-tr{A{x, e~^x)D'^u^) - e~^b{x, e~^x) ■ Du^ + = f{x, e~'^x) in , 

Du^ ■ 7 = g{x, e~^x) on dOe ■ 

As far as the coefficients are concerned, we assume that 

A, b, /, 7 and g are bounded, Lipschitz continuous, Z^-periodic with 

respect to the fast variable and there exists u > such that A > vid. 

In order to formulate the result, we introduce first the associated cell problem which is 
used for the equation inside the domain. We assume that 

'for all x,p E M^, there exists a bounded, Z^-periodic in y, solution 

v{p,x,y) of 

- tr{A{x, y)DlyV) + h{x, y) ■ [p + Dyv) = in 
^ which depends smoothly on {p,x,y). 



(4.i; 



(HI) 



(H2) 



The cell problem for the second corrector is 

for all x, p e and M e , there exists a unique constant -F(M, p, x) 
such that the equation 
(H3) \ -tr{A{x, y) (M + 2{DlyV + MDl^v) + Dl^w) - h{x, y)- 

{Dlv + MDlv + Dyw) = /(x, y) - F{M, p, x) in 
has a bounded solution w, which depends smoothly on {M,p,x,y). 

We remark that in this linear context the assumed uniqueness yields that F is an affine 
function of p and M. 

As it was mentioned in the Introduction, (H2) and (H3) contain two assumptions. One 
is the fact that the cell problems have, for all x, p and M, bounded solutions. The second 
is, of course, the smoothness of v and w with respect to all the variables. This is a 
technical assumption, which is difficult to avoid and, most probably, is true for smoother 
coefficients. 
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We turn next to the boundary condition. The results of Sections |2] and [3] yield, for 
each p,x E M.^ , the existence of a unique constant fi{p, x) such that the boundary ergodic 
problem 

r -tr(A(x, y)DlyZ) - b{x, y) ■ {p + Dyz) = in O , 
(4.2) ) yy 

y {Dz +p) ■-f{x,y) = g{x,y) - i2{p,x) on dO , 

has a bounded, periodic in y', solution z. 
Now we are in position to state our result. 

Theorem 4.1. Assume (HI), (H2), (H3) and that p, satisfies (LI) and (L2). The family 
{u'^)s>o converges locally uniformly, as e ^ Q, to the unique solution u of 




F{D^u,Du,x) + u = in {xn > 0} , 
jl{Du,x) = on {xn = 0} . 



The uniqueness of fj,{p, x) yields that /i is an affine function of p and, hence, of the form 

fi{x,p) = 7(x) -p-gix) , 

for some 7 G C{dO,R^) and g G C{dO). 

The stability and the uniform C*^'"-estimates also yield that the boundary ergodic 
problem has a solution. To prove, however, that 7 and g are Lipschitz continuous is more 
difficult, since it is necessary to take into account the dependence an x of A, b and /, i.e., 
the equation inside the domain. This is why we assume that /2(p, x) satisfies (LI) and 
(L2) in order to have a uniqueness result for (14.31) . 

Proof. Since the proofs that the m^'s converge to viscosity subsolutions and supersolutions 
of (14. 3 p are similar, here we present the argument only for the subsolution case. 
To this end, we introduce the half-relaxed limit 

u{x) = limsup*M^(x) = limsup m^(x') 

and assume that z is a strict maximum of m — </>, where is a smooth test-function. 

If a; G {xiv > 0}, the conclusion follows easily using the perturbed test-function. Indeed 
consider maximum points Xg of 

u^(x) — [(pix) + ev{D(j){x), X, e^^x) + e'^w{D'^(j){x), D(j){x), x, e~^x)) . 

The perturbed test-function is smooth and the proof follows as in the formal asymptotic 
expansion. It is worth pointing out that technically the smoothness of v is required because 
of its dependence on D(f){x) and x is important in the £~^-term. On the contrary, for w, 
the terms are less sensitive and the dependence on x, D(f){x) and D'^(f){x) is enough. 

This is exactly the difficulty we face, when a; G {xat = 0}. Assuming that the solution 
z of (14.21) is smooth with respect to all variables is too restrictive and, most probably, 
uncheckable. We avoid this in the following way. First we write 



z{y] D(t){x),x) = v{x; D(j){x),y) + z{y, x) 
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It follows that z solves the problem 

-tr(A(x, y)Bly£) - b{x, y)-DyZ = inO. 
Dz ■ 7 = ^(x, y) — Ji{D(j){x), x) on dO , 



(4.4) 
where 



9{x, y) = g{x, y) - D(f){x) ■ -f{x, y) - Dyv{x, D(t){x), y) ■ 7(0;, y) 



Since z is not a smooth function of x and, perhaps, we use Theorem 12. II to solve the 
boundary value problem 



I max {-ti{A{x ,y)DlyZ^) -h{x' ,y) ■ DyZ^) = in C ^ 
(4.5) ' 



max (^Dz^ ■ 7 — g{x', y)) = — /i*^ on dO 

\x'~x\<5 



Using the available estimates on z and ft , it is then easy to prove that, as 5 — > 0, 
z^ converges locally uniformly to a solution of (14. 4p and, more importantly, using the 
uniqueness of fi{D(l){x),x), p,^ that converges to fi{D(f){x),x). 

Next we apply the perturbed test-function method in the following way. We look at 
maximum points Xe of 

M^(x) — 0(x) — ev{x, D(f){x),e~^x) — z^{e~^x) — e^w{x, D(f){x), D"^ (t){x) , x) . 

To conclude, it suffices to remark that, since — x as e — 0, we have |xe — a;| < 5 for e 
small enough and (14.51) provides the right inequalities to conclude, since the maximum is 
bigger than the value at x = x^. The fact that z^ may be a non-smooth viscosity solution 
of (14. 5 p creates no real difficulty in the argument. Indeed, it suffices to double variables 
using the test-function associated to (14. 5p . □ 



We remark that v carries all the information to build the second corrector w. On the 
other hand, z appears to be necessary to treat the boundary condition. For this reason, 
our decomposition of the first corrector appears natural. 

We turn next to the nonlinear problem. We use the same notations and, in particular, 
the same domain as in the previous section, and, to simplify the presentation, we consider 
a model problem of the form 




F{eD'^u'^, Du'^,x,e ^x) + e-u^ = ef{x,e ^x) in , 
L{Du^,x,e-^x) = ondOs, 



and assume that F and L satisfy (FO), (Fl), (F2), (F3), (LI) and (L3). Note that, 
compared to the linear setting here, we have multiplied the equation by e here to simplify 
the presentation. 
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Our assumption on the first cell problem is that 

for all x,p E M^, there exists a bounded, periodic in y solution v{p, x, y) 
of 



(H4) 



F{D^yV,p + DyV,x,y) = in 



^ which is a smooth function with respect to all its variable. 

Next we assume that F is smooth in M and p. Using the regularity of v in (H4) we set 
Aip,x,y) = FMiDl v,p + DyV,x,y) and b{p, x,y) = Fp{Dl v,p + DyV, x,y) . 



(H5) 



The cell problem for the second corrector is: 

For all X, p E and M G 5*^, there exists a unique constant 
F{M,p,x) such that the equation 

-tT{A{x, p, y) (M + 2{DlyV + MDl^v) + Dlyw) 



-b{x, p, y) ■ {Dlv + MDp + Dyw) = f{x, y) - F(M, p, x) in 



has a bounded, periodic in y solution w, which depends smoothly on 
M, p, X and y. 

Finally, to specify the boundary condition, we assume that (El) and (F6) are also 
satisfied, and we use the results of Section [2] and [3], which yield the existence of a unique 
constant /2(p, x) such that the boundary ergodic problem 



(4.7) 



F{DlyV,p + DyV,x,y) 
L{p + DyV,x,y) 



in O 

-fM{p,x) 



on dO . 



has a bounded, periodic in y', solution v. 
The result is 



Theorem 4.2. Assume (FO), (Fl), (F2), (F3), (F6), (El), (LI), (L3), (H4), (H5) and 
that F is a smooth function of M and p. If ft satisfies (LI) and (L2), the sequence {u^)e>o 
converges, locally uniformly, as e 0, to the unique solution u of 

F{D'^u, Du,x) + u = m {sAr>0}, 

Ji{Du,x) = on {xn = 0} . 

We omit the proof since it follows very closely the proof of Theorem 14.11 We remark 
that the nonlinear boundary condition does not introduce any additional difficulties. The 
most restrictive assumption concerns the smoothness of F in M and p, which is not easily 
satisfied by Hamilton- Jacobi-Bellman equations. At least Theorem 14.21 applies to cases 
where F is of the form F = Fi + H, where Fi a linear second order operator and H is 
nonlinear in p and smooth in p and x. 
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5. Open Problems and Remarks 

The periodicity of the equation and the boundary is a key assumption in our ap- 
proach. The existence of the boundary cost and of the associated solution of the Neumann 
problem uses this property in a very strong way. 

It would be interesting to know whether approximate correctors exists when the equa- 
tion and the boundary condition are periodic but the half-space domain does not have 
the right periodicity property, i.e., it is typically a hyperplane with an irrational slope. 
This seems to lead to a framework rather close to the one considered by Ishii [13] , but at 
the moment we are unable to obtain results in this direction. Of course the ultimate goal 
is to consider the homogenization problem in random environments. 

The following remark illustrates one of the difficulties to obtain such type of result. For 
q G \ {0}, consider the half-space 

Hq = {x eM.^ : g . a; > 0} 

and assume that the problem 

r -Au = in ifg , 

1^ Du ■ n = g{x) + /i on dHq , 

has a bounded solution u for a continuous, Z^-periodic. Then it is easy to see that we 
must have 

= - lim \B{0,R)nHg\-'^ g{x)dx. 

JB{0,R)nHg 

Now, in M^, choose g{xi,X2) = g{x2) where ^ is a Z-periodic function. It follows 
/^(e2) = .^(0), while, if g = 62 + aei with a is a small parameter, /i(e2) = g{s)ds. This 
shows that a priori fi{q) is not a continuous function of q. Therefore an argument by 
approximation of a non-periodic situation by a periodic one cannot be used to prove the 
existence of fi and of the approximate corrector in the non-periodic framework. 
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